Cluster Dynamics of Planetary Waves 
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The dynamics of nonlinear atmospheric planetary waves is determined by a small number of 
independent wave clusters consisting of a few connected resonant triads. We classified the difi'erent 
types of connections between neighboring triads that determine the general dynamics of a cluster. 
Each connection type corresponds to substantially different scenarios of energy flux among the 
modes. The general approach can be applied directly to various mesoscopic systems with 3-mode 
interactions, encountered in hydrodynamics, astronomy, plasma physics, chemistry, medicine, etc. 



1. Introduction. Planetary-scale motions in the ocean 
and atmosphere are due to the shape and rotation of the 
Earth, and play a crucial role in the problems of weather 
and climate prediction Oceanic planetary waves af- 
fect the general large-scale ocean circulation, can inten- 
sify the ocean currents such as the Gulf Stream, as well as 
push them off their usual course. For example, a plane- 
tary wave can push the Kuroshio Current northwards and 
affect the weather in North America Q]. Atmospheric 
planetary waves detach the masses of cold, or warm, air 
that become cyclones and anticyclones and are responsi- 
ble for day-to-day weather patterns at mid-latitudes [1] . 
Recently a novel model of intra-seasonal oscillations in 
the Earth atmosphere has been developed [3| in terms 
of isolated resonant triads of planetary waves (for wave 
numbers m,£ < 21). The complete cluster structure de- 
ends on the spectral-domain size, both for atmospheric 
and oceanic Q planetary waves. In particular, an 
enlargement of at least some of the clusters is possible, 
with growing of the spectral domain. In some cases this 
yields the energy flux between previously isolated clus- 
ters. To justify the basic model [4|] one has to understand 
whether or not existing clusters are capable to adopt ex- 
ternal energy via this mechanism. In [7[ the isomorphism 
(one-to-one correspondence) between clusters of similar 
structure and corresponding dynamical systems has been 
established. This allows for the study of the dynamical 
behavior of similar clusters. It is shown both in numeri- 
cal simulations 0, Q and in laboratory experiments [loj 
that the dynamics of mesoscopic wave systems does not 
obey statistical description (wave kinetic equations). It 
rather needs a special investigation. 

In this Letter we show that the general dynamics of 
big clusters in mesoscopic systems with 3-modc inter- 
actions is determined by the connection types between 
neighbor triads. In particularly, we analyzed clusters 
of atmospheric planetary waves in the spectral domain 
m,i < 1000 which allowed us to justify the model sug- 
gested in d]. 

2. Triad dynamics. Consider three planetary 
(Rossby) waves with frequencies wi, 0^2 and W3, which 



satisfy the conditions of time and space synchronism: 



(1) 



Ull + UJ2 ^ OJ3, 

mi 4- 7712 = ma, 

\il-h\<h<h+i2, 

< ij, j = 1,2,3, 
7^ ij, i ^.7, hJ = 1,2,3, 
,4+4 + ^3 is odd. 



and ijj ~ m/[£(^+l)]. First three equations correspond to 
three- wave resonance on a sphere while two last equations 
provide non-zero coupling coefficient in the correspond- 
ing dynamical system (see [3l for more details). This 
is the simplest possible cluster that is described by the 
dynamical system 

]3\ = Z B^^Bz^ -^2 ^ ^-^1^3: -^.3 ^ — ZBiB2. (2) 

Here Bj = dBj / dt, Bi ~ aiAi is a time derivative with ai 
being explicit functions of the longitudinal wave numbers 
ij, Aj are modes amplitudes, and Z is the interaction co- 
efhcient which is also some function of the wave numbers. 
The equations ([2]) are symmetric with respect to the ex- 
change of two low-frequency modes 1 <J=^ 2. The mode 
with the highest frequency (which in this paper will be 
denoted by the subscript " 3 ") is a special mode. The 
system ^ has two independent conservation laws 



|B3|2 = {ENi- H)N23/NiN2N3 , 
12 



= 

/i3 = + |B3|2 = {EN2~ H)Ni3/NiN2N3 , (3) 
I12 = hs — I23 = l-BiP 



which are linear combinations of the energy E and en- 
strophy H, defined as by 

E = Ei+E2+E3, NiEi + N2E2 + N3E3 . (4) 

Here Ej is the energy of the j-modc and Nj — £j{£j + 1)- 
The solutions of Eqs. ^ are Jacobian elliptic functions, 
and whether or not its dynamics is periodic is determined 
by the energy in the W3-mode (for details see Q)- 

To understand the dynamics of the energy flow within 
a cluster, the first step would be to initiate a small 
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amount of chosen modes and to study afterwards the en- 
ergy exchange within a cluster. Thus, we begin with dis- 
cussing the evolution of the triad amplitudes with special 
initial conditions, when only one mode is substantially 
excited. If Bi{t = 0) > S2(t = 0) and Bi(t = 0) > 
Bait = 0), then hsit = 0) > hait = 0). The integrals 
of motion are independent of time, therefore /13 3> I23 
at each instant of time and hence \Bi{t)\'^ ^ \B2{t)\'^. 
Moreover, \Bi{t)\'^ |i?3(t)p at every instant. Indeed, 
the assumption \Bi{t)\'^ < |i33(t)P yields /13 ~ I23, 
which is not the case. This means that the wi-mode, be- 
ing the only substantially exited at t = can not share 
its energy with the two other modes in a triad. The same 
is true for the ti;2-mode. In this context we call the modes 
with frequencies toi < uj^ and uj2 < uj^ passive modes, or 
P-modes. 

The conservation laws ([3]) cannot restrict the growth of 
the P-modes from initial conditions when only a;3-mode 
is exited. In this case the P-mode amplitudes will grow 
exponentially: \Bi{t)\ , \B2it)\ (x exp[\ZB3{t = 0)\t] un- 
til all the modes will have comparable magnitudes of the 
amplitudes. Therefore we call the Lo^-mode an active 
mode, or A-mode. The A-mode, being initially excited, 
is capable to share its energy with two P-modes within a 
triad. 

3. Connection types within a cluster. An arbi- 
trary cluster in our wave system is a set of connected 
triads. A cluster consisting of two triads that are con- 
nected via one common mode is called a butterfly; a clus- 
ter of three triads with one common mode is called a 
triple-star. The general dynamics of a cluster depends 
on the type of the connecting mode, which is common 
for the neighbor triads. Correspondingly, wc can distin- 
guish three types of butterflies (with PP-, AP- and AA- 
connections) , four triple-stars (with PPP-, PPA-, PAA- 
and AAA-connections), etc. We begin with considering 
the butterfly and the triple-star dynamics; we then dis- 
cuss an actual dynamics of the more involved but con- 
crete topology of connected triads in atmospheric plane- 
tary waves. 

A PP-butterfly consists of two triads a and 6, with 
wave amplitudes Bja, Bjh, j = 1,2,3, that are con- 
nected via one common mode, Bia ~ Bn,, , which is 
passive in both triads. The equations of motion for this 
system read 

Bla = ZaB^aB^a + ZhB^i^B^b, Bia = Bib, , 
B2a = ZaBl^Bsa , -B26 = ZhB^^B^b , (5) 
Bsa ~ —ZaBiaB2a , B^b = — ZbBiaB2b ■ 

An examination of Eqs. (O shows that they have three 
integrals of motion: 



^23a 

lab ~ 



-\B2a? 



\B: 



3a I 



2b 



\B: 



3a I 



\B: 



(6) 



The first two, 723a and /23b, do not involve the common 
mode Bla = Bib, and are similar to integral /23, Eq- (IS]), 
for an isolated triad. Similarly to the case of the evolu- 
tion of a triad from the initial conditions with an excited 
passive mode, the following conclusion can be made. If 
at t = the amplitudes of one triad substantially ex- 
ceed two remaining amplitudes of the butterfly, that is, 
if I^iqI, |B2a|, |-B3a| > |B2a | , | -83^ | , then this relation per- 
sists. In other words, in a PP-butterfly any of two triads, 
a or 6, having initially very small amplitudes, will be un- 
able to adopt energy from the second triad during its 
nonlinear evolution. 

An AP-butterfly consists of two triads a and b, 
with wave amplitudes that are connected via the common 
mode B^a = Bib, which is active in one triad (a for 
concreteness) and is passive in the second triad (b). In 
this case equations and integrals of motion arc: 



- ZaB2aBza , B^b — —ZbB3aB2b , 

- ^a-t>lQi>3a , i>2b — ^b^3a^3b , 
: —ZaBiaB2a + ZbB^bB^b , 

Il2a = \Bla\'^ — \B2a\'^ , hsb = |-B2b|^ + IB^- 



33b 



In 



\Bla\ 



\B. 



\B 



(7) 



(8) 



3b 



The integrals /12a and /23b, do not involve a common 
mode B^a Bib] they are similar to the corresponding 
integrals I12 and 123, for the isolated triad. In the case, 
when the triad a is excited at i = much stronger then 
the 6-triad (in which case /12a ^ ^236) then the smallness 
of the positively definite integral of motion /23b prevents 
the triad b from adopting energy from the triad a during 
all the evolution. The situation is different, when triad 
b is excited initially and /23b ^ Ii2a- In this case the 
initial energy of the triad b can be easily shared with the 
triad a. The smallness of /12a only required that during 
the evolution \Bia\ ~ |i?2a|- 

For an AA-butterfly with a common active mode in 
both triads {B^a — B^b) we have: 



Bla ~ ZaB2aB3a , Bib = —ZbB^^bB^a 
B2a = ZaBl^Bsa , -S2b 
-S3a = —ZaBiaB2a — ZbBibB2b 



ZbBl^B^a , (9) 



Il2a — |7?laP 
lab = l-BlaP -\ 



-\B2aV 
\B3a? 



h 



2b 



\Bib?-\B 



2b 



\B 



(10) 



3b 



3b| 



Again, the integrals /i2q and /12b do not involve a com- 
mon mode B^a ~ Bib and are similar to I12, Eqs. ([3]), 
for an isolated triad. Simple analysis of these integrals of 
motion shows that the energy that is initially hold in one 
of the triads will be dynamically shared between both 
triads. 

Finally we consider one of the triple-star clusters, e.g. 
the APP-star, in which a common mode is active in the a- 
triad and passive in the b- and c-triad: B^a = Bib = Bi^. 
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Clust. 


M 


Modes [m,£] 


Af2 


Connecting triads 


Ai 


1 


[4,12] [5,14] [9,13] 




- 


A2 


2 


[3,14] [1,20] [4,15] 


2.1 


[4,15] [10,24] [14,20] 








2.2 


[1,20] [14,29] [15,28] 








2.3 


[1,20] [15,75] [16,56] 


A3 


3 


[6,18] [7,20] [13,19] 


3.1 


[2,15] [5,24] [7,20] 


A4 


4 


[1,14] [11,21] [12,20] 


4.1 


[1,14] [9,27] [10,24] 




5 


[2,6] [3,8] [5,7] 


5.1 


[4,14] [9,27] [13,20] 




6 


[2,6] [4,14] [6,9] 








7 


[6,14] [2,20] [8,15] 


7.1 


[2,20] [11,44] [13,35] 


CXl7,8 


8 


[3,6] [6,14] [9,9] 


7.2 


[2,20] [30,75] [32,56] 








7.3 


[32,56] [26,114] [58,69] 


1x19,10 


9 


[3,10] [5,21] [8,14] 


— 


— 




10 


[8,11] [5,21] [13,13] 








11 


[2,14] [17,20] [19,19] 


11.1 


[2,14] [18,27] [20,24] 




12 


[1,6] [2,14] [3,9] 


11.2 


[6,44] [14,21] [20,24] 




13 


[3,9] [8,20] [11,14] 


11.3 


[9,35] [11,20] [20,24] 


Kl 11-16 


14 


[1,6] [11,20] [12,15] 


11.4 


[3,20] [45,75] [48,56] 




15 


[9,14] [3,20] [12,15] 








16 


[2,7] [11,20] [13,14] 







TABLE I: In the first 3 columns the following data in the do- 
main m,^ < 21 are given: cluster's form, triad numbers and 
modes within a cluster; in the last two - numbers of connect- 
ing triads and their modes are given that enlarge correspond- 
ing cluster when spectral domain m,i < 1000 is regarded. 



This system has four integrals of motion: 

Slap — |-S2aP , 



Il2a 

l23b = \B2b\^ + \B3b\\ 
labc = l-BlaP + I^SaP 



'23c 



IB: 



2c 



>3c 



(11) 



\B: 



36 1 



3c I 



Similarly to butterflies, there exists one integral of mo- 
tion for each connected triad, that does not involve the 
common mode: these are the integrals /12a, hsb, and 
^23c, which are the same as the corresponding integrals 
in the isolated triad. The integrals /23b and /23c prevent 
the b- and c-triads (which are connected via a P-modc) 
from adopting energy from the initially excited a-triad. 
In those cases when the b- and/or c-triads are initially 
exited, the a-triad can freely adopt their energy via the 
connecting A-mode. 

Any triad that is connected to a cluster (no matter 
how big the cluster would be) via its passive mode cannot 
adopt energy from the cluster, if the triad is not excited 
initially. On the other hand, a triad that is connected to 
any cluster via an active mode can freely adopt energy 
from the cluster during its nonlinear evolution. 
4. Topology and cluster dynamics for atmospheric 
planetary waves. The frequencies of atmospheric plan- 
etary waves arc 




11.4^ 




14^S/^12 


\13 / 


\ 16 / 


|11^3^><ri1.1 

/ 11. 


\ 11 / 



^2nmj/e.j{i, + 1) . 



(12) 



FIG. 1: Color online. Triads belonging to butterfly 1x12 (left) 
and to caterpillar (right) are drawn by bold (red) lines while 
new, connected to them triads appearing in spectral domain 
m,l < 1000 are drawn by thin lines. (Red) arrows are coming 
from an active mode and show directions of the energy flux. 
The numbers inside each triangle correspond to the numera- 
tion in Table 1. 



The negative sign indicates wave propagation opposite to 
the rotation of the Earth (east to west), opposite to the 
Earth rotation (with frequency O). The integers €j with 
and rrij < £j describe the eigen-mode structure, which 
is j-spherical harmonics with £j and £j — rtij zeros in 
the longitudinal and latitudinal directions, correspond- 
ingly. The Diophantine Eqs. ()12p and lui+uj2 = t^^a, have 
many solutions, each of them describing an exact reso- 
nance of ideal planetary waves (ignoring the real Earth 
topography, etc.). In this approximation we can describe 
all important resonant triads in the Earth atmosphere, 
see Table 1. First we restricted oursclf to the so-called 
meteorologically significant spectral domain with wave 
numbers £,m < 21 [J. In this domain we have found 
in 0] four isolated triads, denoted as A1...A4, three 
PP-butterflies txi5fi, ^7.8i and ixig.io, involving of six tri- 
ads A5 . . . Aio and one, further called caterpillar, Klu-ig, 
consisting of six triads An . . . Aig with three PP-, one 
AP- and one AA-connection. 

In this Letter we show that the topological structure 
of clusters in the extended spectral domain m,£ < 1000 
is richer, in particular, some clusters are enlarged by new 
resonances formed by modes with m,£ > 21. To illus- 
trate this mechanism, we introduced following rows in 
the Table I: 1) JVi - the number of clusters in the spec- 
tral domain m,£ < 21, with the same numeration as given 
m [3; 2) "Modes" - resonance clusters belonging in the 
same spectral domain as in 1); 3) J\f2 - number of the 
additional clusters in the spectral domain m,£ < 75; 4) 
" Connecting triads" - additional resonance clusters which 
appeared in the spectral domain as in 3). The topologi- 
cal structure of resonance clusters in the spectral domain 
TO, £ < 75 is shown in FigUl where numbers inside of tri- 
angles correspond to the numeration A/i , 7V2 in the Table 
1. 

There exist altogether 1965 isolated triads and 424 
clusters consisting of 2 to 3691 connected triads, among 
them 235 butterflies, 95 triple-triad clusters, etc. - see 
the histogram in Fig. ^ The three largest clusters con- 
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2 4 6 8 10 12 14 16 

Size of cluster 



FIG. 2: Color online. Horizontal axes denotes the number of 
triads in the cluster while vertical axes shows the number of 
corresponding clusters. 

sist of 14, 16 and 3691 connected triads. For the clarity 
of presentation we did not display on the histogram the 
largest 3691-cluster, which we will call further the mon- 
ster. 

It can be seen that 82.2% of all clusters are isolated 
triads. Their dynamics has been investigated in [J in 
all details. The energy oscillates between three modes 
in the triads, the period of this oscillation being much 
larger than the wave period. It was found to be inversely 
proportional to the root-mean-square of the wave ampli- 
tude. 

235 clusters ( ~ 10.5%) are butterflies. Among them 
there are 131 PP-, 69 AP- and 35 AA-buttcrflies. The 
butterfly dynamics is restricted by three integrals of mo- 
tion, and it can be shown that the phase space of butter- 
flies is four-dimensional. In 0] only PP-butterflies have 
been considered. For the initial conditions studied in this 
Letter their dynamics is similar to that of two isolated 
triads. 

Preliminary numerical simulations show that only in 
the case when the initial levels of excitation in both triads 
are compatible, a periodic energy exchange is observed 
not only within a triad but also between two connected 
triads. 

The 95 triple-triad clusters include 66 linear clusters 
with two pair connections, 25 "3-stars" with triple con- 
nections and 4 -triangles (with three pair connections). 
A similar classifications can be performed for all the re- 
maining clusters. For example, the monster includes one 
mode (218,545), participating in 10 triads, three modes, 
participating in 9 triads, 5 modes - in 8 triads, 23 - in 7, 
50 - in 6, 90 - in 5, 236 ~ in 4, 550 - in 3, and 1428 modes 
- in 2 triads (butterflies). The analysis of their dynami- 
cal behavior depends critically on the type of connection, 
as was shown above. 

For example, the 16-triad cluster can be divided 
into "almost separated" parts (connected through PP- 
connections) parts: 5 triads, one AA- and one AP- 
butterfly, one AAP-star and one AAA-star with an AP- 
conncctcd triad. The overall qualitative conclusion is 



that even big clusters are dynamically not very different 
from separated small clusters with active connections, 
AA-butterflies, AAA-stars, etc. 

To clarify the dynamics of the first sixteen triads 
A1...A16, it is important to establish their (possi- 
ble) connection to the clusters in a bigger spectral do- 
main. We have found that the triad Ai remains iso- 
lated, A3 turned into an isolated PA-butterfly. It can 
be proven [ll| that these objects remain "forever" iso- 
lated, even if the size of the bigger domains goes to in- 
finity. The triads A2 and A4 became parts of the mon- 
ster but are connected with it via P-modes; in this sense 
they are practically separated. Three initial butterflies, 
IXI5 g, [XI7 8, ocig 10 have PP-connections and thus their dy- 
namics does not differ much from the dynamics of an 
isolated triad. Moreover, an increase of the spectral do- 
main to m,^ < 1000 does not change the situation sub- 
stantially: the butterfly ixig 10 remains isolated in an ar- 
bitrarily large (even infinite) domain, IXI7 g became part 
of a 5-triads cluster, also with PP- and PPP-type of con- 
nections, (see Figdl), 1x15 is now part of the monster but 
only via PP-connections. The caterpillar gains one P- 
connected triad (11.4 in the Tab. 1) and one AAA-star 
(11.1, 11.2 and 11.3 triads) with two P-connections, see 
FigUl Therefore almost all 16-triads, (except of the AA- 
butterfly ixiis^ig) can be considered as completely or as 
almost separated from the rest of the system. 
Conclusions: 

• In the physically relevant domain of atmospheric 
planetary waves (m , £ < 1000, when the mode-scale is 
larger than the height of the Earth atmosphere) we have 
determined and described topology of all clusters that are 
formed by resonantly interacting planetary modes. The 
cluster set contains isolated triads and sets of 2-, 3-, . . . , 
16 and 3691 connected triads, with 2-3-, . . . , 9-mode 
and (maximum) 10-mode connections. 

• Analyzing the integrals of motion we suggested a 
classification i) of triad modes into two types - active 
(A) and passive (P), and ii) of connection types between 
triads - AA, AP and PP. We have shown that in AA- 
butterflies the energy can flow in both directions, in AP- 
butterflies only from one triad to the second one (and 
not vice versa), while in PP-butterflics the triads are "al- 
most isolated" . We have also shown that the dynamical 
behavior of bigger clusters can be similarly characterized 
by connection types like AAA, AAP,..., etc. 

• As a first approximation, almost all triads in the me- 
teorologically significant domain of Tabl. 1, can be con- 
sidered as completely or almost separated from the rest of 
atmospheric planetary waves, and therefore energy oscil- 
lations within them can lead to intra-seasonal oscillations 
in the Earth's atmosphere, as suggested in 

• Our analysis is based on the structure of the dy- 
namical system ([2]) and has two advantages: i) it does 
not need to exploit explicit form of the interaction co- 
efficients Z; ii) it is completely analytical without using 
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numerics (which will be required for more general ini- 
tial conditions). Thus our results can be used directly 
for arbitrary resonant 3- wave systems governed by these 
equations, e.g. drift waves, gravity-capillary waves, etc. 

• A touchstone of any theory is, of course, an exper- 
iment. Numerical experiments with resonance clusters 
described by barotropic vorticity equation are now on 
the way, preliminary results confirm our theoretical con- 
clusions. The next step of utmost importance would be 
to study some physical mechanisms that might destroy 
clusters. Speaking mathematically, one can always intro- 
duce big enough resonance width 

17 = 01 + 02 - 03 > 

such that a substantial part of cluster's energy will be 
re-distributed among other waves via non-resonant inter- 
actions (see Fig. 2, [12|). From the physical point of view, 
the source of the resonance broadening might have sub- 
stantially different reasons - from baroclinic instabilities 
due to the effects of the free surface at large scales to the 
effects due to the Earth topography at smaller scales to 
the increasing the level of turbulence, say, in summer due 
to increasing sun activity, and thus going into the regime 
of fully developed wave turbulence, to the inclusion of 
dissipation and forcing. This effects can also be com- 
bined, of course. The problem of the utmost importance 
is therefore to study the resonance clusters behavior in 
the situation when at least some of these effects are in- 
cluded. The analysis applied to climatic variations on 
geological scales "typically give indications of low dimen- 
sionality and (...) the hope of justifying the modelling 
of weather or climate in terms of a small set of ordinary 
differential equations" [l^l. But it does not mean, of 
course, that the overall energy flux will stay nicely reg- 
ular. Indeed, as it was shown in [l^, a special choice 
of instabilities included into ([2|) will cause appearance 
of strange attractors. The study of this transition from 
regular to chaotic regimes in resonant clusters and mu- 
tual interrelations of relevant physical parameters is the 
subject of our further research. 

• Last not least. Whereas it is quite difficult to check 
experimentally the theory for planetary waves, it can be 
done much easier within a framework of laboratory ex- 
periments with some water waves, e.g. gravity-capillary 
waves. Some preliminary program of laboratory exper- 
iments with this type of waves has been worked out in 
[3] but only for simplest triad's clusters for at that time 
the algorithms of cluster computing [l^, [3l were not 
developed yet.) Presently this program can easily be 
elaborated for clusters of some more complicated struc- 
ture. Considerably more sophisticated preliminary work 
is needed for planning laboratory experiments with grav- 
ity water waves. As it was shown in [l^l, a 3- wave res- 
onance system differs principally from wave systems in 
which 4- and more wave resonances are allowed. Indeed, 
in a 3-wave resonances system, each nonlinear resonance 



generates new scale, while already in a 4- wave system 
this is not necessarily true. Indeed, beginning with 4- 
wave resonances, different types of energy fluxes have 
been pointed out: scale-resonances (as in 3-wave system), 
angle-resonances (formed by two couples of wavevectors 
with pairwise equal lengths) and mixed cascades. In this 
case, clusters have more complicated structures presented 
explicitly in [l8l | for 4-wave interactions of gravity water 
waves. Qualitative dynamics of small clusters formed by 
resonant quartets is briefly as follows. "One wave mode 
may typically participate in many angle resonances and 
only one scale resonance. Thus, one can split large clus- 
ters into "reservoirs", each formed by a large number of 
angle quartets in quasi-thermal equilibrium, and which 
are connected with each other by sparse links formed 
by scale quartets." (see [3|). This means, in particular, 
that scale-resonances cause spectrum anisotropy. On the 
other hand, one can easily compute what ratio aspect of 
the sides of a laboratory tank should be chosen in order 
to suppress scale-resonances. In this case, regular pat- 
terns on the water surface are to be expected, similar to 
what was observed in BSHH. It would be interest- 
ing to study these experimental data in order to establish 
these nearly permanent patterns observed, which can be 
attributed to some specific resonance clusters. Existence 
of independent resonance clusters can shed some light 
on the origin of Benjamin- Feir instabil ity or McLean 
instability [2^ (see also Discussion in [fsjl. 

Another interesting and even more complicated area 
of further investigations would be the study of the dy- 
namical behavior of resonant quintets, sextets and so on. 
Methods developed in 1^, l3| allow to compute the clus- 
ters and corresponding wave frequencies. This informa- 
tion can be afterwards used, for instance, for investigat- 
ing of some special types of shallow water instabilities 
2J|. In systems containing simultaneously quartets and 
quintets [25| , the explicit construction of clusters can pos- 
sibly yield the explanation of the competing regimes be- 
tween these two types of instabilities. We point out here 
three possible scenarios due to the change of 1) bound- 
ary conditions, 2) frequency range under the study, and 
3) initial distribution of energies among the modes of 
a cluster. In the first two cases, some quartets can be 
suppressed, thus turning the quintets into the principal 
clusters (and vice versa, of course). In the third case, 
a lot of different sub-scenarios are possible: quartet and 
quintet connected via one common mode can become " in- 
dependent" if energy in the common mode is too small; 
a quartet can " die out" if its initial energy is not enough 
for nonlinear interactions; same for the quintet; energy 
can flux from the quartet to the quintet via the common 
mode, in dependence on the energy distribution among 
the other modes in the quartet thus yielding the dying 
out of the quartet; same for the quintet; etc. Combin- 
ing the methods [l^, [lB| and results of [i^ will help to 
single out some of the possibilities immediately. For in- 
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stance, one can compute wavenumbers correspondin g to 
a resonant cluster, put them into the Eq.(2.11) from |25 | 
and check whether coupling cociJicicnt(s) in this cluster 
is (are) non-zero. Another advantage would be construc- 
tion of the topological structure of the complete cluster 
set as it is done in [l^ instead of resonance curves: it 
gives a general overview of the nonlinear dynamics rather 
than a (possible) subset of resonances formed with a fixed 
wavevector (in particular cases this subset can be empty, 
of course). One more advantage would be the following. 
As it was shown in [2^, instability analysis is based on 
the properties of the roots of the polynomial given by 
Eq.(3.10) (see also Fig.3a-3d therein). In some cases a 
qualitative instability analysis can help but "to get qual- 
itative results in general case one should solve (3.10)" 
([2^, p. 313). This computations are very involved while 
the polynomial has degree 4 and its coefficients depend on 
the dynamical coupling coefficients. On the other hand, 
knowledge of explicit wavenumbers for modes forming a 
resonance cluster will turn the coefficients of (3.10) into 
known constants and the problem can be easily solved. 
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